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Several theories have been proposed to describe the formation of black hole seeds in 
the early Universe and to explain the emergence of very massive black holes observed 
in the first thousand million years after the Big Bang’ *. Models consider different 
seeding and accretion scenarios* ’, which require the detection and characterization 
of black holes in the first few hundred million years after the Big Bang to be validated. 


Here we present an extensive analysis of the JWST-NIRSpec spectrum of GN-z11, 
an exceptionally luminous galaxy at z= 10.6, revealing the detection of the 
[NeIv]A2423 and CII*A1335 transitions (typical of active galactic nuclei), as well 

as semi-forbidden nebular lines tracing gas densities higher than 10° cm”, typical 
of the broad line region of active galactic nuclei. These spectral features indicate 
that GN-z11 hosts an accreting black hole. The spectrum also reveals a deep 

and blueshifted CIVA1549 absorption trough, tracing an outflow with velocity 
800-1,000 km s”, probably driven by the active galactic nucleus. Assuming local 
virial relations, we derive a black hole mass of log(Mgy/Mo) = 6.2 + 0.3, accreting 
at about five times the Eddington rate. These properties are consistent with both 
heavy seeds scenarios and scenarios considering intermediate and light seeds 
experiencing episodic super-Eddington phases. Our finding explains the high 
luminosity of GN-z11 and can also provide an explanation for its exceptionally high 


nitrogen abundance. 


GN-z11 was recently observed with JWST. The analysis of the NIRCam 
images revealed an unresolved nuclear component and a disc-like 
component with a few 100 parsec (pc) radius®. A first NIRSpec spec- 
trum was presented in ref. 9, which found it to be consistent with star 
formation, although the presence of an active galactic nucleus (AGN) 
was not excluded. Here we explore the latter scenario using a deeper 
spectrum of GN-z11. 

Figure 1a shows the detection of the [NeIv]AA2422,2424 doublet. 
As NelV requires photons more energetic than 63.5 eV, this line is an 
unambiguous AGN tracer”? and not seen in star-forming galaxies, 
not even those hosting the Wolf-Rayet (WR) stars”. 


We also detect CII*A1335 emission (Fig. le). This line is commonly 
observed in AGN**. In star-forming galaxies, this line is generally 
totally undetected; when detected, it is extremely faint and always 
associated with deep CIIA1334 resonant absorption”, not seen in GN-z11. 

The [NIv]A1483, NIv]A1486 doublet is very sensitive to the gas density 
(and insensitive to ionization parameter, metallicity and shape of the 
ionizing spectrum). Figure 1c shows the detection of the semi-forbidden 
NIVJA1486 line (critical density 4.7 x 10° cm”) and the non-detection 
of the forbidden [NIV]A1483 line (critical density 1.5 x 10° cm™), which 
indicate densities much higher than 10° cm”. Specifically, the various 
photoionization models shown in Fig. 2a (see the Methods for details) 
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Fig. 1| Magnification of the spectra of GN-z11 around specific spectral 
features of interest, along with their single and multiple Gaussian models 
(Methods). Dashed lines indicate the rest-frame wavelengths of the lines at 
z=10.603. a, [NeIV]AA2422,2424 doublet; b, NIII] multiplet, illustrating the 
detection of the resolved N111]A1754 emission; c, NIV] doublet, showing the 
absence of [NIV]A1483 despite the strong NIV]A1486; d, CIV blueshifted 
absorption trough and redshifted resonant emission, compared with the CIV 


illustrate that the upper limit on the doublet ratio requires densities 
greater than about 10° cm”, which are incompatible with the densities 
of the ionized interstellar medium (ISM) that are typically in the range 
of 10-10? cm”, and only rarely approach a few times 10* cm” (ref. 18). 

Even stronger constraints come from the NIII] multiplet (Fig. 1b). 
This is contributed primarily by four semi-forbidden lines at 1,748.6 A, 
1,749.7 A, 1,752.2 A and 1,754.0 A (and the much weaker 1,746.8 A). 
All of these have high critical densities (>10° cm), but the 1,748.6 A 
and 1,754.0 A transitions have the highest critical density of 10° cm™ 
(note that atomic physics requires fixed flux ratios F,754/F,743 = 1.05 and 
Fya6/ F759 = 0.14; Methods). The 1,754.0 A line is well resolved compared 
with the rest of the multiplet, and its intensity is well constrained to be 
0.25 + 0.04 of the total intensity of the multiplet. This high ratio can be 
achieved only when both the 1,749.7 A and 1,752.2 A transitions are sup- 
pressed relative to the 1,748.6 A and 1,754.0 A because of the very high 
density. As shown in Fig. 2b, when compared with the expectations from 
awide range of photoionization models, the observed ratio requires a 
density higher than 10° cm” at 30 (higher than 10° cm” at 20). These 
high densities are completely inconsistent with any HII regions in any 
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P-Cygni profile observed in low-metallicity, young star-forming galaxies 
(stack: orange dashed line; most extreme case: orange dotted line), showing 
inconsistency with the latter. e, CII/CII*AA1334,1335 doublet (seen in emission, 
without P-Cygni, only in type 1 AGN); f, Expected flux of the NIV1718 line in the 
case that NIV]1486 was associated with WR stars. In a-c, e andf, the continuum 
is subtracted, whereas in d the continuum is normalized to one. The grey 
dotted lines indicate the noise level (10). 


star-forming galaxy but are fully in the realm of the broad line regions 
(BLRs) of AGN, which are characterized by extremely high densities 
(about 10°-10" cm”). 

Therefore, the most plausible explanation is that GN-z11 hosts an 
AGN and that these semi-forbidden lines observed in its spectrum are 
mostly emitted by the associated BLR. 

It may seem puzzling that the spectrum of GN-z11 does not seem 
to show the typical ‘broad lines’ seen in type 1 quasars and AGN, with 
widths of thousands km s”. However, the width of the broad lines scales 
quadratically with the black hole mass, hence in the case of ‘small’ black 
holes the broadening is expected to be substantially smaller. Moreover, 
there are classes of type 1 AGN that have broad lines with widths less 
than1,000 kms": these are the so-called Narrow Line Seyfert 1 (NLSy1), 
the permitted lines of which are broader than their forbidden lines, but 
not by a large factor, and in many cases reaching a width of only a few 
100 km s” (ref. 19), and which are inferred to have small black holes 
(about 10°M,) (ref. 20). This seems to be the case of GN-z11, in which the 
semi-forbidden (NIII], NIV]) and permitted lines (Mgll) all have widths 
between 430 kms ‘and 470 km s”, whereas [Nelll] has a significantly 
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Fig. 2 | Flux ratios of density-sensitive nitrogen lines as a function of 
hydrogen gas density, n,,. A large range of Cloudy models (Methods) are 
compared with the values observed in GN-z11. Models with metal-poor 

(Zeb = 0.1Z,) and metal-rich (Z,,., =1Z,) gas are shown with solid lines and 
dashed lines, respectively (colour-coded according to the ionization parameter 
U), inthe scenario in which either an AGN (filled symbols demarcating different 
black body temperatures for the accretion disc, Taon) or stellar populations 
(open markers for various ages, t.) is responsible for the incident radiation 
field. a, [NIV]A1483/NIV]A1486 flux ratio. b, Ratio of NIIIJA1754 to total flux of 
the multiplet. The black dashed lines and blue-shaded regions (in decreasing 
darkness for 10, 20 and 30 confidence level as indicated) show the observed 
fractional contribution of N1I1JA1754 and upper limit on [NIV]A1483/NIV]A1486 
obtained for GN-z11, indicating that the gas emitting these lines has high 
density (n, 2 10° cm at 30). The light-green-shaded areas highlight the range 
of densities typical of the broad line regions (BLRs), whereas the grey-shaded 
regions highlight the range of densities typical of the ionized interstellar 
medium (ISM). 


narrower width (340 + 30 kms‘), hence coming from the host galaxy 
(either from HII regions or from the narrow line region of the AGN). We 
note that the interpretation of some permitted and semi-forbidden 
lines, such as the Balmer lines and CIII] is made complex by the fact 
that these are also generally contributed to by the ISM photoionized 
by star formation in the host galaxy. 

The spectrum of GN-z11 also shows a deep (EW,. = 5 A) and 
blueshifted absorption trough of the CIVAA1548,1550 doublet 
(Fig. 1d). Deep CIV absorption is sometimes observed in young stel- 
lar populations, but the depth observed in GN-z11 would require high 
metallicities, typically solar or super-solar”". This is in contrast with 
the metallicity inferred from the nebular lines of GN-z11 (Z = 0.1Z,) 
(ref. 9). To illustrate more quantitatively the inconsistency with the 
stellar-wind origin, the orange dashed line in Fig. 1d shows the stacked 
spectrum of local galaxies with metallicity around the value inferred 
for GN-z11, resampled to the NIRSpec grating resolution: the stel- 
lar trough is much shallower than that observed in GN-z11 and with 
acompletely different shape. Apart from the stellar origin, this deep 
CIV absorption is seen also in lower redshift star-forming galaxies and 


associated with galactic outflows”. However, in these cases the outflow 
velocities are only of afew 100 km s” (Methods), whereas for GN-z11 the 
CIV trough traces a much faster outflow of around 800-1,000 km s”. 
A more plausible explanation of the deep blueshifted trough of CIV 
is that GN-z11 is part of the class of broad absorption line (BAL) qua- 
sars, which are characterized by deep absorption of blueshifted CIV 
by up to several thousands km s”. Actually, GN-z11 would fit in the 
‘mini-BAL’ category, with velocities between 500 and 2,000 kms, 
more common in lower luminosity AGN, or in the ‘narrow (approxi- 
mately 1,000 km s”) absorption line (NAL) quasars category”. The 
spectrum also shows aclear CIV redshifted emission, whichis probably 
tracing the receding component of the outflow. As CIV is a resonant 
line, this is the counterpart of the redshifted Lya identified in ref. 9 
(consistent shift and width). 

In sum, the detection of [NelV] and CII*, the extremely high gas 
density matching those of the AGN BLRs, and the presence of a deep, 
blueshifted absorption trough of CIV tracing a high-velocity outflow 
are all consistent with the scenario in which GN-z11 hosts an accreting 
black hole, thatis, an AGN, specifically what would be called NLSyl and 
(mini-)BAL/NAL AGN. 

Inthe Methods, we also discuss other diagnostics, suchas the ratio of 
UV transitions (for example, CIII]/CIV, CIII]/Hell) and the upper limits 
on high ionization lines (NVAA1238,1242 and [Nev]A3426), are fully 
consistent with the AGN scenario. 

Some works have suggested that GN-z11 may host a population of WR 
stars”. The HellA1640 line shows a potentially broad profile, as shown 
in Extended Data Fig. 1b (although the wings are mostly in the noise). 
This, if confirmed, could come from the inner region of the BLR but 
could also trace the presence of a WR population. However, various 
other features are inconsistent with the main contribution from WR 
stars. Specifically, inthe case of WR stars, the NIVAA1483,1486 doublet, 
if present, is always accompanied by an even stronger NIVA1718 line, 
with a prominent P-Cygni profile, which is not seen at high confidence 
in the spectrum of GN-z11 (ref. 13) (Fig. 1f); [NelV] and CII* are never 
seen associated with WR stars”; when present, the NIII] multiplet has 
a much weaker 11754 component”. Therefore, if WR stars are present 
in GN-z11, then they must co-exist with the AGN and are unlikely to 
play a dominant part in the excitation of the observed nebular lines. 

Assuming local virial relations, the black hole mass can be esti- 
mated fromthe line widths and continuum luminosity. As discussed in 
the Methods, we estimate a black hole mass of about 1.6 x 10°M,. 
In the Methods, we also discuss potential uncertainties and caveats 
in the determination of the black hole. 

We infer a bolometric luminosity of the AGN of 10* erg s” (Methods), 
which is a factor of about 5 higher than the Eddington limit (with an 
uncertainty of a factor of 2). Super-Eddington accretion is generally 
inferred for NLSyls and is one of the scenarios proposed to rapidly 
grow supermassive black holes in the early Universe*”. 

Figure 3 shows how the black hole mass in GN-z11 would have evolved 
at earlier cosmic epochs if accreting at the Eddington rate, or at the 
super-Eddington rate estimated at the time of observation. For com- 
parison, the grey-shaded areas show the range of possible black hole 
seeds scenarios: black holes resulting from the direct collapse of pri- 
mordial clouds into seeds with masses in the range of about 10*-10°M,, 
the so-called direct collapse black holes (DCBHs); rapid merging of stars 
and black holes in dense, nuclear star clusters; accretion onto Popula- 
tion III black hole seeds or even normal stellar remnants!*°. Many of 
these semi-analytical models and cosmological simulations could 
reproduce the mass of GN-z11 at z= 10.6 (refs. 7,25-27). The solid and 
dashed lines show the evolutionary tracks for some of them (described 
more extensively in the Methods). These can be broadly divided in 
models assuming heavy seeds (DCBH), whose accretion is limited to the 
Eddington rate, and intermediate mass (stellar clusters) or light (stellar 
remnants) seeds experiencing episodes of super-Eddington accretion. 
It is interesting to also note that GN-z11 evolving at sub-Eddington 
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Fig. 3 | Black hole mass as a function of redshift (on a logarithmic scale) and 
age of the Universe. The black hole mass inferred for GN-z11 is shown with the 
large golden symbol. The red-shaded region indicates the evolution expected 
in the case of super-Eddington accretion at the level inferred for GN-z11. The 
darker-blue-shaded region shows the black hole mass evolution assuming 
Eddington-limited accretion, whereas the lighter-blue-shaded region shows 
the case of evolution in the case of sub-Eddington accretion (between 0.1land1 
the Eddington rate). The horizontal grey-shaded regions indicate the range 

of black hole seeds expected by different scenarios. Solid and dashed lines 
indicate the evolutionary tracks of various simulations and models””** that can 
reproduce the GN-z11 black hole mass, with different seeding and accretion 
rate assumptions, as detailed in the Methods. The small grey symbols indicate 
the black holes measured in quasars (QSOs) at z = 6-7.5 (refs. 1,2) (whose 
representative loerror bar is shown inthe top left), most of which can originate 
froma progenitor suchas the black hole in GN-z11. 


rate can easily result into the supermassive black holes (107-10°M,) 
observed in quasars at z= 6-7, as predicted by many models. 

However, it is possible that the local or low-z scaling relations do 
not apply for AGN at such early epochs. If we disregard the local virial 
relations and instead assume that the black hole in GN-z11 is accreting 
at the Eddington rate, then the black hole mass would be 10’M,. A black 
hole with this mass is more difficult to account for, but achievable 
in models assuming heavy seeds and episodes of super-Eddington 
accretion’”°”’, 

Taking the stellar mass of the extended disk-like component meas- 
ured in ref. 8 (M, = 8 x 108M,,), it is possible to locate GN-z11 on the 
Moy-Mgrar relation. As shown in Fig. 4, GN-z11 is placed above the local 
relation, although marginally consistent within the scatter. An early 
evolution above the local Mg,-M,4, relation is what is expected from 
models invoking DCBHs and/or super-Eddington accretion?*”. 

We note that the exceptionally high nitrogen abundance inferred for 
GN-z11 (specifically, high N/O)? becomes much less problematic inthe 
AGNscenario. To begin with, several ‘nitrogen-loud’ AGN have already 
been found both at low and high redshifts (including NLSy1)”**". So 
GN-z1l is not a very peculiar system in this context. Second, the mass 
of the BLR in AGN is very small*: 


Luy n = 
Moir = Lela erg alor ) Mo 


Specifically, the Hy luminosity observed in GN-z11 (L,, = 1.7 x 10” erg s*, 
even assuming the extreme case it is not contributed by the NLR and HII 
regions) implies a mass of the BLR of only a few solar masses. It would 
take just one or two supernovae to enrich such a small mass to solar 
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Fig. 4 | Black hole versus stellar mass diagram. The location of GN-z11 (large 
goldensymbol) is compared with local galaxies as indicated by the small red 
symbols and their best-fit relation (black solid line and uncertainty traced by 
the grey-shaded region)*°. The grey symbols show the values estimated for 
quasars (QSOs) at z = 6-7 (ref. 37), although in these cases the galaxy mass is 
inferred from dynamical tracers. The blue symbols are AGN at z > 4 for which 
the black hole and galaxy stellar mass has been measured withJWST data 
(Methods) using the same calibration as in ref. 36 for consistency. 


or super-solar metallicity, especially within the small physical region 
associated with the BLR (about 10° pc for GN-z11) (ref. 33). These could 
be supernovae from supermassive stellar progenitors, with high nitro- 
gen yields**. However, even without invoking exotic scenarios, given 
the accelerated metal enrichment of sucha small, central region, it is 
also possible that secondary, recycled nitrogen production can occur 
within a timescale ofa few tens Myr (especially given the very fast cool- 
ing times at such high densities, which allow star formation to quickly 
occur out of cooled SN ejecta). 

We finally discuss our results within the context of the recent JWST 
findings of an excess of exceptionally luminous galaxies at high redshift. 
GN-z11is one of the first of such hyperluminous galaxies at high-z to 
be spectroscopically confirmed, and for which such a detailed spec- 
troscopic analysis has been feasible. The AGN scenario revealed by our 
analysis provides a natural explanation for the exceptional luminosity 
of GN-z11. If this is representative of the broader class of luminous gal- 
axies discovered at high-z, then it would greatly alleviate the tension 
with models and simulations. 
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Methods 


Observations and data processing 

The data presented in this paper are part of the JADES survey* and, 
specifically, obtained through programme ID 1181 (principal investi- 
gator D. Eisenstein). GN-z11 was observed in two epochs: the first one 
on UT 5 and 7 February 2023 and the second one on UT 4 and 5 May 
2023. The February observations were already presented in ref. 9. We 
refer to that paper for a detailed description. Briefly, the spectroscopic 
data were obtained with four different configurations of the NIRSpec 
micro-shutter array (MSA)? ™, using athree-shutter nodding pattern. 
Four different dispersers were used to cover the 0.6-5.3 um wavelength 
range: the low-resolution prism mode (exposure time of 6,200 s, per 
configuration), and three medium-resolution gratings (3,100 s each, 
per configuration), which provide a nominal spectral resolution of 
R =1,000 fora uniformly illuminated slit”. However, the highly compact 
light profile of GN-z11, with respect to the width of the slit, results in 
a substantially higher effective resolution. To estimate the effective 
resolution, we forward model the morphology of GN-z11 through the 
NIRSpec instrument for the grating dispersers, finding that the resolu- 
tion ranges between 1,100 and 2,100. Four MSA configurations were 
used (two pointing and two dither positions). The May observations 
were similar, but in this case they consisted of three consecutive dith- 
ers (with three different MSA configurations), each with three nods, 
resulting in an on-source exposure of 2.7 h for the prism, each of the 
three medium-resolution gratings and also with the high-resolution 
grating G39SH/290LP. Unfortunately, at the location of GN-z11 onthe 
MSA the latter spectrum is heavily truncated at wavelengths longer 
than 4.1 um, hence all strong optical emission lines are not observed 
with this grating. 

By combining the two sets of observations, the total exposure time 
is 9.6 h with the prism and 6.15 h with each of the medium-resolution 
gratings. 

The data processing is also described in ref. 9, and we refer to that 
paper for a detailed discussion. Here we only mention that we used the 
pipeline developed by the ESA NIRSpec Science Operations Team and 
the NIRSpec GTO Team. Most of the processing steps in the pipeline 
adopt the same algorithms used in the JWST Science Calibration Pipe- 
line”. Different from the official pipeline, the final one-dimensional 
(1D) combined spectra are obtained by combining the 1D individual 
spectra rather than performing the extraction process in the combined 
two-dimensional spectra. This step guarantees that the final 1D spectra 
are well flux calibrated for slit losses. In the combination process, we 
also applied a 30-clipping algorithm and excluded bad pixels based 
onthe data quality files provided by the pipeline. The extraction of 1D 
spectra in the individual exposures is also optimized on the basis of 
science. In this paper, we adopt a three-pixels (0.3”) extraction along 
the slit, as it improves the signal-to-noise ratio (S/N) of the spectrum 
(for point sources). Finally, the GTO pipeline provides spectra beyond 
the nominal wavelength range for the spectral configuration G140M/ 
FO7OLP by taking into account the transmission filter throughput in 
the flux calibration processing step. The extended spectra cover the 
wavelength range of 1.27-1.84 um. 

Moreover, here we combine the grating spectra in their overlapping 
ranges, which increases the S/N in those regions. The combined spectra 
in these regions were resampled to 8 A around the NIV doublet, not to 
affect resolution, and to 12 A around the CIV, as in this case higher S/N 
is required on the continuum to properly trace the CIV absorption. 

Inthe paper, we adopt the flat ACDM cosmology from Planck18 with 
H = 67.4 km s” Mpc ‘and Q,, = 0.315 (ref. 43). 


Emission line fitting 

The emission lines were fitted using single or multiple Gaussian lines 
and a simple power law for continuum subtraction. The best-fit para- 
meters for the continuum and Gaussian components were found using 


the MCMC (Markov chain Monte Carlo) algorithm to estimate the 
uncertainties. For the purposes of this paper, each line in the rest-frame 
UV was fitted independently, except for doublets or multiplets, whose 
line widths were forced to the same value (but see discussion below 
for the CIII] doublet) and the relative wavelength separation of the 
doublet and multiplet was forced to the nominal rest-frame wave- 
length. The absolute velocity of each line (or group of lines in the case 
of doublets and multiplets) was not constrained to the exact redshift 
given in ref. 9, to allow for small wavelength calibration uncertain- 
ties associated with the positional uncertainties of the target within 
the shutter. (The uncertainties in the target acquisition may result in 
the target being offset by up to about 0.05” relative to the nominal 
position, which is used by the pipeline for the wavelength solution. 
This unknown offset would introduce a wavelength offset by up to 
0.5 spectral pixel, which corresponds to a different velocity offset 
(given the wavelength-dependent resolution) in different regions of 
the spectrum.) (Jakobsen, personal communication) and also to allow 
for small velocity shifts between different lines, which are commonin 
AGN, and especially in the BLR. We restricted our fitting to the lines of 
interest for this paper. 

Extended Data Table 1 provides alist of the fitted emission line widths 
and fluxes and Extended Data Fig. 1 shows the additional fitted lines 
not shown in the main text. 

In the case of the NIII] multiplet, the 1,748.6 A and 1,754.0 A transi- 
tions come from the same upper level, hence their flux ratio is fixed 
by the associated Einstein coefficients, specifically Fy754/F,743 = 1.05. 
Similarly, 1,746.8 A and 1,752.2 A come from the same upper level and 
their flux ratio is fixed to Fy746/F,75. = 0.14. The inferred line widths are 
deconvolved from the line spread function as inferred for the GN-z11 
light profile. 

For the CIII]1906,1908 doublet, it is not possible to resolve the two 
components; attempting to fit it with two components makes the fit 
degenerate between the width and intensity of the two components. 
The additional caveat of this CIII] doublet is that it is also commonly 
seen in normal star-forming galaxies and in the NLR of AGN, so it can 
also have a contribution from the host galaxy, as for the [Nelll] emis- 
sion. In ref. 44, the authors use IFS spectroscopy to reveal that the CIII] 
emission is resolved on scales of several 100 pc. As a consequence, we 
do not include the (spectrally) unresolved CIII] in our analysis, as it 
does not provide constraints on either the BLR or the host galaxy. In 
Extended Data Table 1, we report the total flux and width using a single 
Gaussian. However, for the sake of completeness, we report that by 
fitting two components with separate full-width at half maximum 
(FWHM), accounting for the NLR and BLR, gives C1I1]A1906/A1908 of 
0.627129 for the narrow components, with FWHMs of 314 +120 kms 
(consistent with the [Nelll] width), and 560+ 80 km s” for the 
C11JA1908 broad component (consistent with the NIV] width). 

The [011]3726,3729 doublet would potentially be an additional for- 
bidden line, detected in the observed wavelength range, which could 
be used to constrain the velocity dispersion in the host galaxy. However, 
unfortunately, the doublet is unresolved. Attempting to fit it (by forc- 
ing the two components to have the same width) results ina FWHM of 
365 +55 kms ‘anda flux ratio of 0.627031. 


CIV absorption and emission 

In this section, we provide some additional details on the CIV absorp- 
tion. As mentioned in the text, CIV P-Cygni profiles with a significant 
CIV blueshifted trough are seen associated with atmospheres of young, 
hot stars. Yet, the depth of this feature is a strong function of metalli- 
city”, and the deep trough observed in GN-z11 would require stars with 
solar or evensuper-solar metallicities, completely inconsistent with the 
much lower metallicity inferred for GN-z11. To illustrate the inconsis- 
tency with the pure stellar origin, we have stacked 11 UV spectra from 
the CLASSY HST survey”, with metallicity around the value inferred for 
GN-zllinref. 9 (Z=0.1Z,), specifically 7.4 < 12 + log(O/H) < 7.9. We were 


conservative by excluding galaxies with strong CIV emission. We also 
excluded one WR galaxy, as we discuss that the spectrum cannot be 
dominated by WR stars (see main text and section ‘The WR scenario’). 
The continuum of the spectra was normalized to one by using a simple 
linear fit in the spectral ranges 1,410-1,480 A and1,560-1,600 A, con- 
sistent with the analysis of the spectrum of GN-z11 in the same spectral 
region (Fig. 1d). The resulting stacked spectrum is shown witha dashed, 
orange line in Fig. 1d and illustrates inconsistency with the trough seen 
in GN-z11. To be conservative, in Fig. 1d, we also show the case of the 
most extreme spectrum among the 11 selected, the one with the deep- 
est CIV absorption. Although the wings of the stellar winds can extend 
out to 2,000 km s7, the profile and depth observed trough at these 
metallicities is inconsistent with the observed trough in GN-z11 at So. 

The blueshifted CIV trough (and redshifted emission), therefore, 
is not a P-Cygni feature associated with stellar (atmosphere) winds. 
Rather, it is tracing a galactic outflow, as observed in lower redshift 
starbursts” and in (mini-) BAL/NAL AGN’. The determination of the 
velocity requires knowledge of the exact wavelength of the redshifted, 
rest-frame CIV transition. Unfortunately, there are small wavelength 
uncertainties associated with each grating because of the uncertainties 
of the location of the sources within the shutter, as discussed above. In 
this specific case, we calibrate the velocity shift based on the NIV line, 
which is in the same gratings and has a similar ionization potential as 
CIV. The outflow velocity is also subject to different definitions. The cen- 
troid of the trough relative to the mean of the two CIV transitions gives 
avelocity of -790 km s”. If we consider the blue edge of the trough rela- 
tive to the bluest of the two transitions (CIVA1548.19), then we obtaina 
velocity of -1,040 km s”. These velocities are significantly higher than 
those inferred from the CIV absorption in starburst-driven outflows””°, 
but in the range of BAL quasars that can span from 500 kms “to several 
thousands km s” (refs. 23,45-49,51,52). 

The classification boundary between mini-BAL and NAL AGN is not 
sharp, with different authors giving different definitions in terms of 
width and/or blueshift of the absorption”*** ”. Here we simply give 
a generic classification as mini-BAL/NAL without aiming at a more 
specific category. 

It should be noted that some past works have reported some rare 
starburst galaxies showing outflows with high velocities, evenin excess 
of 1,000 km s” (refs. 53-56). However, these outflows are traced by 
lower ionization transitions (Mgll absorption and [OIl] emission). More 
importantly, a close inspection of those cases reveals that each of them 
shows some AGN signature ([NeV] emission and/or broad Mgll emission 
and/or broad HB emission and/or X-ray emission and/or located inthe 
AGN or composite region of diagnostic diagrams). Therefore, although 
the AGN contribution to the bolometric luminosity of these galaxies 
may be arguable (also taking into account the variable nature of AGN), 
itis likely that the high-velocity outflows seen in these rare cases are 
actually driven by the AGN that they host. 

Finally, it should be noted that the CIV absorption trough goes nearly 
to zero (asin many BAL quasars), which implies the total covering factor 
of the emitting source by the outflowing ionized gas along our line of 
sight. However, the errors leave scope for a contribution of 30% of the 
emission potentially not covered by the CIV absorption, which can be 
associated with the extended host galaxy. Yet, ifhigher S/N data confirm 
the CIV trough going to zero, this would imply that the outflow has an 
extent covering also the host galaxy, that is, about 400 pc, which would 
be fully consistent with recent findings of BAL outflows extending on 
scales of up to several kpc (refs. 49,57-59). 

Given that CIV is a resonant line, the observed redshifted emission 
is also tracing the CIV counterpart of the redshifted Lya emission seen 
in ref. 9—that is, the receding side of the outflow. 

We finally note that the spectrum between Lya and the NV doublet 
shows the tentative signature of an NV blueshifted trough (Extended 
Data Fig. 1d), which would be associated with the highly ionized outflow, 
but it requires additional data to be confirmed. 


Constraints from other emission lines and diagnostics 

Although the paper focuses ona few lines discussed in the main text, in 
this section we also discuss other emission lines that have either lower 
S/N, more severe blending or whose upper (or lower) limits give line 
ratios that are fully consistent with the AGN scenario. 


Møll and CIII]. The MglI2796,2804 doublet is well resolved with the 
grating and in principle a good tracer of gas density in the range be- 
tween 10° cm and 10“ cm”. However, the observed ratio, 1.3670:%%, is 
so uncertain to be consistent both with the low-density regime (ratio 
of around 1) and the high-density regime (ratio of about 2). Moreover, 
even if additional data allows constraining the MglI doublet ratio more 
tightly, these are resonant transitions, which are, therefore, strongly 
sensitive to the optical depth and radiative transfer effects. 

The C111JAA1907,1909 doublet would also be a good density tracer, as 
the ratio of its two components is primarily sensitive to the gas density 
and changes strongly between 10* cm and 10° cm? (with the blue com- 
ponent A1907 going to zero at high densities), similar to the NIV] doublet. 
However, as discussed above, the two components are unresolved with 
the grating, and we cannot obtain reliable constraints on the gas density 
or onthe line widths. More importantly, CII] emission is commonly seen 
also in star-forming galaxies and in the NLR of AGN, so it may partially 
come also from the low-density ISM of the host galaxy, as is the case for 
[Nelll]. As already mentioned, recent IFS observations show CIII] to be 
resolved on scales of several 100 pc (ref. 44). Itis interesting that when 
fitted with narrow and broad components, as discussed in the previous 
section, the narrow component gives widths formally consistent withthe 
[Nelll], whereas the broad componentis consistent with the NIV width. 


NV and NeV. Additional transitions from species requiring ionizing 
photon energy higher than about 60 eV, suchas NV and NeV (in addition 
to NelV seen in GN-z11), are often seen as evidence for the presence of 
an AGN. Yet, conversely, their absence should not be necessarily seen 
as evidence for the absence of an AGN, as often these lines are weak 
even in AGN and remain undetected if the S/N is not high enough” ®. 
Moreover, the intensity of these lines varies strongly from case to case. 

With the prism it is not possible to assess the presence of NV because 
it is blended with Lya and its damping wing. Regarding the gratings, 
the G140M band, in which NVis redshifted, is the least sensitive of the 
three medium-resolution spectra. Although there is a hint of the NV 
doublet (Extended Data Fig. la 20 integrated signal) we obviously do 
not quote it as a tentative detection. The inferred upper limit onthe NV 
emission is not very constraining, but the important aspect in the con- 
text of this paper is that itis still fully consistent with the presence of an 
AGN. We demonstrate this in Extended Data Fig. 2, in which the upper 
limits on the NV/CIV and NV/Hell ratios for GN-z11 are compared with 
a sample of the broad lines in type 1 AGN and also with a sample of 
the NLR in type 2 AGN™, and illustrating that the non-detection of NV 
is fully consistent with the AGN scenario. 

It is also interesting to compare NV with NIV, as this ratio is not 
dependent on the nitrogen abundance, although NIV is detected (or 
reported) less frequently in AGN. In the well-studied type 1.8 AGN at 
z=5.5, GS-3073 (refs. 16,30,64), the NV is five times fainter than NIV, 
which would be totally undetected in our spectrum. Inthe type 1 quasars 
exploredin ref. 65, the NIV broad line is very strong, whereas NV is unde- 
tected, with an upper limit that is about 10 times lower than the NIV flux. 

NeVis also not detected, neither in the grating nor in the prism spec- 
trum. The upper limit on the NeV/Nelll ratio is about 0.2. However, in 
ref. 63, the authors have shown that AGN models can have NeV/Nelll as 
low as 107-10 *. Hence the non-detection of NeV is also not constrain- 
ing about the presence of an AGN. 

Finally, we note that AGN accreting at super-Eddington rates havea 
lower energy cutoff, and hence are less likely to emit hard photons that 
can produce highly ionized species, such as NV and NeV. 
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Hell and CIV. Hell is detected in the prism and, more marginally, inthe 
grating (Extended Data Fig. 1). 

As already discussed, CIV is detected in the grating, but with a P-Cygni 
profile, hence its flux is a lower limit because of self-absorption. 

The interpretation of these limits using photoionization models 
is very much model-dependent. We illustrate this in Extended Data 
Fig. 3a,b. Specifically, Extended Data Fig. 3a, as in ref. 9, shows the 
location of GN-z11 on the CIII]/CIV versus Hell/CIII] diagram and in 
which the red-squared and blue-starred symbols show the location 
of models from refs. 10,66 for the NLR of AGNs and for star-forming 
galaxies, respectively, and ina range of about +0.3 (see legend) dex of 
the metallicity inferred in ref. 9 for GN-z11. GN-z11 can be consistent 
with both AGN and star-forming models. 

Extended Data Fig. 3b shows the same diagram in which we instead 
plot the models from ref. 67, inthe same (low) metallicity range for both 
AGN and SF galaxies. In this case, GN-z11 is much more consistent with 
the AGN models and inconsistent with the models for star-forming 
galaxies. 

Yet, if the permitted and semi-forbidden lines of GN-z11 are coming 
from the BLR, as argued in this paper, then neither of the models above 
actually apply, as they are developed for the low-density environments 
of the NLR and HII regions. It is, therefore, more instructive to com- 
pare with the line ratios observed in the BLR of type 1 AGN. These are 
taken from the compilation of ref. 61 and shown with purple circles 
in Extended Data Fig. 3c. The line ratios observed in GN-z11 are fully 
consistent with the broad lines of type 1 AGN. For completeness, inthe 
same panel we also plot the ratios observed for the NLR of type 2 AGN, 
compiled in ref. 62 (mostly overlapping with the ratios observed 
for the broad lines), and the star-forming galaxies from the CLASSY 
survey". 


The WR scenario 
In this section, we discuss the scenario recently proposed that GN-z11 
may be similar to local WR galaxies”. 

The Hell marginal detection shows a potentially broad profile (about 
10,00 km s”, although the broad wings are mostly in the noise), which 
may be associated with the inner BLR, but also may resemble the broad 
Hell profile characteristic of WR stars. Therefore, there might be a 
contribution from WR stars and, specifically, WN stars, given the strong 
nitrogen lines. 

However, there are various spectral features that cannot be 
accounted for in the WN scenario. 

WN stars are also characterized by very strong NIVA1718 resonant 
emission, stronger than the NIVA1486, and typically with a prominent 
P-Cygni profile”. In GN-z11, despite the very strong NIVA1486, there is 
no trace of the NIVA1718 line. Figure 1f shows the spectrum of GN-z11 at 
the expected location of NIVA1718 and in which the shaded red region 
shows the expected intensity of the line, based on the strength of the 
NIVA1486 line. The GN-z11 spectrum is totally inconsistent with the 
presence of the NIVAI718 WR signature. 

Furthermore, neither [NeIv]A2424 nor CII*A1335 are ever seen associ- 
ated with the WR population”. 

Finally, even if WN show prominent NIII] emission, the strength of 
the A1754 component of the multiplet is much fainter in WR galaxies 
suchas Mrk966 (ref. 17) and consistent with densities typical of the ISM. 

In sum, although WR stars might be present in GN-z11, they are 
unlikely to dominate the excitation of most nebular lines. 

Extended Data Table 2 summarizes more schematically the observa- 
tional features consistent or inconsistent with the AGN scenario, the 
WR scenario and a compact starburst without WR stars. 


Photoionization modelling 

We used the CLoupy photoionization code® to explore the effect of 
varying physical conditions on some emission line ratios constrained 
byJWST/NIRSpec. The primary goal is to explore the ratios of emission 


lines within a given doublet or multiplet, hence lines of the same ion 
(specifically NIII and NIV) that are effectively insensitive to the chemical 
abundance and ionization parameter, while sensitive to density and 
only with secondary dependence on temperature. For this reason, the 
details of the photoionization models are not as important as when 
exploring other line ratios. We considered a nebula of constant pres- 
sure in plane-parallel geometry. However, we have verified that other 
scenarios, suchas a cloud with constant density, do not affect our find- 
ings. For completeness, we considered both AGN and stellar templates 
for the shape of the incident radiation field. Its normalization is set 
by the ionization parameter, defined as U = ®,,/(n,c), where @, is the 
surface flux of hydrogen-ionizing photons at the illuminated face of 
the nebula, ną is the number density of hydrogen and cis the speed of 
light. The hydrogen density and ionization parameter were varied in 
logarithmic steps of 1, respectively from ny =1cm° upton, =10" cm”, 
and starting at log, ,U=-—3 and ending at log, ,U = -1 (refs. 10,66,67). 

In the AGN scenario, we adopted the multi-component continuum 
template implemented in CLOUDY, consisting of a black body anda 
power law, varying the black body temperature (Ton = 10° K and 10° K) 
while fixing the power-law slope to a = -1.4 (note that this is the slope 
underlying the black body at energies above the Ly-edge) and leaving 
other optional parameters as default. For the AGN models, we consid- 
ered gas-phase metallicities of Zae =0.1Z, and Z,,., =1Z,. By contrast, 
the star-forming models are restricted to Ze = 0.1 Zo, as the hard ion- 
izing spectra of metal-poor stars are essential to form sufficient triply 
ionized nitrogen (requiring 47.5 eV), whose presence in GN-z11 is evi- 
denced by the strong NIV emission (EWyyy ag = 9-0 + 1.1 A; ref. 9), whereas 
metal-rich stars would not produce enough hard ionizing photons to 
make the NIV line visible. In the star-formation scenario, we used stellar 
population synthesis models, including binary stars generated by BPASS 
v.2.1 (ref. 69) for a single burst of star formation (with varying ages, 
t,/Myr € {1, 10, 100}), assuming the same metallicity as the gas (that is, 
Z. = Zeb = 0.1Z,) and an IMF” that ranges in stellar mass from 1M, to 
100M,,. Both inthe AGN and star-formation cases, calculations are run 
until a neutral hydrogen column density of N = 10” cm” is reached to 
ensure that in all models the nebula is matter bounded; we note, how- 
ever, that the highly ionized nitrogen lines are produced in the very inner 
part of the cloud, such that the boundary conditions do not significantly 
affect our results. In total, this results ina parameter grid of 15 different 
densities, 3 ionization parameters, 3 temperatures or stellar ages, 2 or 
1 metallicities for the AGN and star-formation models, respectively, 
or a total of 15 x 3 x (3 x 2 + 3 x 1) = 405 possible model configurations. 

The relevant nitrogen line ratios for all of these (except for eight cases 
in which CLOUDY reported a failure) are shown in Fig. 2, from which we 
conclude that they are consistent between the AGN and star-formation 
scenario, and their density dependence is largely independent of ioniza- 
tion parameter, metallicity or the precise shape of the incident radiation 
field (that is, AGN or star formation and the corresponding parameter 
Tion OF 6). 

At high densities, the NIVA1483/NIVA1486 ratio approaches zero 
(ny 2 10° cm), whereas N111]A1754 plateaus at a fractional contribution 
to the multiplet of about 0.23 at higher densities still (n4 z 10 cm’), 
both pointing towards the presence of a broad line region in GN-z11 
given the observed values. 

Finally, to increase the readability of Fig. 2, we have separated the AGN 
and star-forming models in two separate panels in Extended Data Fig. 4. 


Continuum shape 

If GN-z11 is a type 1 AGN, then we should be directly seeing the light 
from the accretion disc. In the case that the accretion disc dominates, 
the UV-to-optical continuum should follow a simple power law of the 
form F, « A’ with B = -7/3 = -2.33 (ref. 71), as observed in type 1 AGN, 
and NLSy1””’, modulo the UV turnover whose wavelength increases 
with black hole mass and also modulo effects of dust reddening, which 
often makes the spectrum redder. 


Inthe case of GN-z11, the spectrum is contributed to also by the under- 
lying galaxy identified in ref. 8 inthe NIRCam images. This compo- 
nent is significantly fainter than the nuclear point-like component. It 
is difficult to quantitatively establish its contribution to the spectrum, 
because part of the light may fall outside the shutter, and ina different 
fraction inthe four dither or pointing positions, and not easy to recon- 
struct because of the slight positional uncertainties discussed above. 
In Extended Data Fig. 5a, we show the contribution from the galactic 
component (dotted-orange line) to the spectrum, assuming that the 
entire light of the galaxy is captured by the spectrum, corresponding 
to about one-third of the flux and using the spectral template inferred 
inref. 8 for the extended component. 

The additional component to take into account is the nebular con- 
tinuum associated with the BLR (as well as any other ionized gas in the 
host galaxy). The BLR typically has alow covering factor”, therefore the 
nebular continuum is not expected to be strong, but its contribution 
must be quantified. In most physical conditions typical of the ionized 
gas in the BLR, NLR or HII regions, the nebular continuum is linked to 
the intensity of the Balmer lines. We have estimated the nebular emis- 
sion using a Cloudy model with a metallicity of 0.1 Z, and a density of 
10° cm” (between the BLR and ISM origin scenarios) and normalized 
to have the same Hy flux as observed in the spectrum of GN-z11. The 
nebular spectrum does not change drastically as a function of density, 
except obviously for the emission of the forbidden and semi-forbidden 
lines; however, our focus is on the nebular continuum, so we ignore the 
mismatch of the emission lines, as a detailed photoionization modelling 
of their flux is beyond the scope of this paper. We note that the nebular 
continuum is also included in the model spectrum fit to the extended 
componentin ref. 8; therefore, not to include it twice, we have measured 
the Hy flux in the ref. 8 spectrum and normalized the Cloudy nebular 
spectrum only to the Hy flux obtained by the difference between the 
observed value and the flux in the ref. 8 model spectrum. The resulting 
nebular spectrum is shown witha dashed purple line in Extended Data 
Fig. 5a. Again, the mismatch of the emission lines should be disregarded, 
as the goal is not to reproduce them with the Cloudy model. 

Extended Data Fig. 5b shows again the observed spectrum, in log-log 
scale, in which the galactic and nebular components have been sub- 
tracted. Although the noise is large, especially at long wavelengths, also 
as a consequence of the model subtraction procedure, the resulting 
spectrum is well fitted by a simple power law, in the parts not affected 
by the emission lines. The best-fitting slope is -2.26 + 0.10, hence con- 
sistent with the continuum expected from an accretion disc. Note that 
this is not evidence in support of the presence of an AGN, as also young 
galaxies may have power-law shapes, it is only meant to show consist- 
ency with the AGN scenario. 

We finally note that, although with a large scatter, the UV spectrum 
of AGN often shows a Fell hump between about 2,300 A and about 
3,100 A (refs. 75-79). The prism spectrum of GN-z11 does not show an 
obvious Fell bump, although a more detailed analysis and modelling 
is required to assess the presence or absence of sucha bump, whichis 
deferred to a separate paper. However, we note that at such early epochs 
thereis little time for the ISM to be enriched with iron through the SNla 
channel®, so a weak or absent Fell bump would not be unexpected. 


Variability 

The luminosity of AGN can be variable, froma few per cent toa factor of 
afew, on short (days) and long (years) timescales. We have investigated 
the possible presence of variability. Before the recent NIRCam images 
obtained in February 2023 (ref. 8), deep photometric observations were 
obtained with HST about 10 years earlier®”, corresponding to about 
l year in the rest frame of GN-z11. Most of the HST photometric data 
points have error bars that are too large to be useful for constraining 
variability. However, the photometric point reported in ref. 81 with the 
F160W filter has a relatively well constrained value: 150 + 10 nJy, within 
an aperture of 0.35”. NIRCam does not have the same filter, however, 


the photometry obtained in the F150W filter can be used and transposed 
tothe F160W filter by using the NIRSpec prism spectrum. After extract- 
ing photometry from a 0.35” aperture (as in ref. 81), and extrapolating 
withthe NIRSpec spectrum, we obtain a F160W equivalent photometry 
of 141 + 2 nJy, which is consistent with the HST previous photometry 
within lo. If we consider that about 30% of the flux is diluted by the host 
galaxy, the comparison of the photometry between the two epochs 
would indicate a variability of 10% at only lo. This is certainly nota 
detection of variability, but it is consistent with the range of variabi- 
lity amplitudes observed in NLSy1 and, more broadly, in type 1 AGN®. 


X-ray emission 

GN-z11 is not detected in X-rays. Flux limits are obtained from the 
Chandra Deep Field North, which was a2 Ms observation performed in 
2002 (see* for final results). Their sensitivity map gives a point source 
limit in the soft (0.5-2 keV), hard (2-7 keV) and full (0.5-7 keV) bands 
of 1.54 x10, 7.9 x 10” erg cm” s and 4.9 x 10” erg cm” s”. Source 
detection requires ano-source probability P < 0.004. The tightest limit 
inthe soft band translates toa rest frame 5.8-23.2 keV luminosity limit 
atz=10.6 of 2.2 x 10” erg s”. Assuming a typical NLS1 photon index of 
2.3 means that Lx (2-10 keV) is less than 3 x 10” erg s”. 

The bolometric correction for NLS1 in the 2-10 keV band, BCy, is 
about 100 (ref. 85). There is a significant systematic uncertainty here 
due to the unseen flux in the FUV, in which the emission is expected to 
peak (see fig. 3 in ref. 86). Moreover, the 2-10 keV flux entirely origi- 
nates from the corona, the early development of which and possible 
dependence on black hole spin are unknown (ref. 87 cautions against 
using his X-ray BC values for NLS1). Proceeding with BC, = 100 means 
that the Chandra upper limit is almost three times above the luminosity 
inferred from the JWST flux at 1,400 A. We predict a conservative SB 
flux of 5 x 108 erg cm™ s”. This would be detectable in about 1 Ms with 
the candidate NASA Probe mission AXIS. The coronal emission from 
local NLSIs is highly variable and the above BC represents a mean value 
(note that the intrinsic disc flux seen in the UV is much less variable®*). 


Black hole mass estimate 

For the vast majority of high redshift AGN, the black hole masses are 
inferred using single-epoch measurements and the so-called virial 
relations, that is, relations between the black hole mass, the width of 
the lines of the BLR and the continuum or line luminosity”’** **. These 
relations are calibrated on nearby AGN, using either reverberation 
mapping techniques and/or direct dynamical measurements of the 
black hole. The black hole mass scales about as the square power of 
the width of the BLR lines and about as the square root power of the 
luminosity, with a proportionality constant that depends on the specific 
waveband (or line) for the luminosity estimation. 

The most accurate virial relations would be those using Ha and HB. 
In our case, Hy could be used as a proxy. However, as discussed, the 
Balmer lines are probably contributed to by the star formation in the 
host galaxy, hence not reliable to trace the black hole mass. 

The CIII] doublet is also sometimes used to infer the black hole mass. 
However, this is not well resolved and, as for the case of the Balmer 
lines, this is probably contaminated by the ISM and star formation in 
the host galaxy. 

Mgll is often used. In our case, the Mgll doublet is clearly detected, 
but the S/N is fairly low for the measurement of the width (Extended 
Data Fig. 1). If we take the width resulting from the fit and the relation 
provided in ref. 95: 
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then we get a black hole mass of 1.4 x 10°M,,. However, given the low 
S/N on the Mgll doublet, we prefer to use as representative width of 
the BLR lines the profile of the high S/N and isolated NIV line. If we 
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adopt this width into the equation above, we obtain a black hole mass 
of 1.6 x 10°M,. The uncertainty is totally dominated by the scatter in 
the virial scaling relation, which is about 0.3 dex (ref. 96). 

Moreover, there are various other systematic uncertainties and cave- 
ats that can affect the black hole mass estimate. To begin with, it is not 
obvious that the local virial relations apply at high redshift. The main 
issue is whether the dependence of the BLR radius on luminosity evolves 
with redshift or not. The most plausible scenario is that the square 
root dependence of the BLR radius from luminosity is primarily set 
by the dust sublimation radius. In ref. 97, the authors argue that, given 
the extremely high densities in the nuclear region of AGN (hence high 
optical thickness even at very low dust-to-gas ratios), unless the nuclear 
region is totally devoid of dust, the same R,,-L relation is unlikely to 
evolve with redshift. Assessing whether the virial relations depend 
on the accretion rate or not is more problematic. On the one hand, 
in ref. 98, the authors argue that the effect of radiation pressure is to 
reduce the effective gravitational force on the clouds of the BLR; the 
net result is that the standard virial relations applied to BHs accreting 
close to the Eddington rate could underestimate the black hole mass by 
a factor of several. On the other hand, reverberation mapping of AGN 
accreting at super-Eddington has revealed that in these cases the size of 
the BLRisa factor of several, and up to an order of magnitude, smaller 
than expected from the R,, -L relation for sub-Eddington AGN (ref. 99 
and references therein), which would imply that the standard virial rela- 
tions overestimate, by a factor of several, the black hole masses in AGN 
accreting at super-Eddington. Overall, it is possible that the radiation 
pressure effect and the offset from the Rgır-L relation might cancel 
each other out. However, currently it is not really possible to provide 
an accurate assessment on how much AGN accreting at or beyond the 
Eddington rate might deviate from the standard virial relations. 

Finally, the black hole masses from other JWST studies at z = 4-8 
(refs. 30,93,97,100,101) are shown in Fig. 4. These are based on the Ha 
or HB width and flux. We clarify that these are re-estimated by using 
the same calibrations used in ref. 36 for local galaxies. 


AGN bolometric luminosity estimate 

We derive the bolometric luminosity of the AGN by using the continuum 
luminosity at A,,., = 1,400 A and the luminosity-dependent bolometric 
correction given in ref. 87: 
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We also assume, as discussed in the previous sections, that 30% of the 
continuum flux at this wavelength is because of the underlying galactic 
component? and that, therefore, the AGN continuum luminosity at 
this wavelength is 0.7 of the observed value. We infer a bolometric 
luminosity of 1.08 x 10* erg s”. The resulting ratio between bolometric 
and Eddington luminosity is 5.5, also affected by an uncertainty of a 
factor of at least 2, coming from the uncertainty on the black hole mass. 


Comparison with cosmological and hydrodynamical 
simulations 
There is a vast literature discussing the formation of early black 
holes and on how they evolve in the first thousand million years, 
by using hydrodynamical and cosmological simulations, as well as 
semi-analytical models. The production and elaboration of models in 
this area have recently seen surge with the goal of specifically interpret- 
ing the results fromJWST. It is beyond the scope of this paper to provide 
an exhaustive description of the assumptions and results of the several 
models and simulations. However, in this section, we briefly discuss 
that many of them can explain the properties of GN-z11 and provide 
some possible constraints on the seeding scenarios. 

We start by considering the results obtained in ref. 7 from the FABLE 
hydrodynamical, cosmological simulation, in which they focused on 


the largest halo at z = 6 (with a virial mass Myo. = 6.9 x 10M, of the 
Millennium box). The latter may appear an extreme choice; however, 
we note that GN-z11 does live in an overdense region and probably at 
the core of a protocluster®’™. In the FABLE simulation, the black hole 
seed has a mass of 10°M, at z= 13. The accretion rate is capped to Edding- 
ton and uses the Bondi—-Hoyle-Littleton-based formalism; however, 
as small scale, non-isotropic accretion is unresolved in the simulation, 
FABLE, like Illustris, uses a Bondi-Hoyle-Littleton rate boosted by a 
factor of 100. Feedback energy in FABLE scales as 10% of the available 
accretion energy, È = eMc?, where e = 0.1is the radiative efficiency of 
the accretion flow. At high redshifts, this is primarily injected as thermal 
energy in the vicinity of the black hole, with a duty cycle of 25 Myr. We 
overplot the fiducial model in ref. 7 in Fig. 3 (orange solid line, labelled 
as B23), illustrating that this can easily reproduce the mass of the black 
hole in GN-z11 at z= 10.6. 

To explain the most massive BHs observed atz = 6-7, the same study 
as above’ also explores the scenario of earlier seeding (z = 18) and allows 
the black hole to accrete at up to two times the Eddington limit; in this 
case, the model could explain a black hole nearly five times more mas- 
sive than GN-z11 at z= 10.6. 

Inref. 35, the authors explored the early evolution of black holes using 
the TRINITY cosmological empirical model!®, which is based on halo 
statistics from N-body simulations and incorporating empirical galactic 
scaling relations. The authors specifically explore the case of GN-z11. 
They illustrate that its mass and black hole to stellar mass ratio can be 
explained by their model starting with an intermediate mass seed of a 
few times 10° seeded at z= 15, accreting on average at sub-Eddington 
rates, but intermittently also at super-Eddington. Their track is shown 
with a solid-teal line in Fig. 3 (labelled as Z23). 

Recently, in ref. 25, the authors have explored the properties of GN-z11 
within the context of the semi-analytical model CAT. They find that the 
black hole mass of GN-z11 and its location on the Mgu-M,ta diagram can 
be interpreted both in terms of light seeds (at z = 20-23) that can have 
super-Eddington accretion phases, or Eddington-limited heavy seeds 
formed at z= 14-16. Out of their various tracks, Fig. 3 shows only two 
samples of their tracks, in the case of a light (red-solid) and a heavy 
seed (red-dashed), which can both reproduce the mass of GN-z11 at 
z=10.6 (labelled as S23). In both cases, the semi-analytical model can 
also reproduce the black hole to stellar mass observed in GN-z11. 

Inref. 3, the authors suggested that the detectability of accreting BHs 
at high redshift by JWST implies that these are probably originating 
from heavy seeds. Specifically, their models can reproduce the mass 
of GN-z11 at z= 10.6 but only with seeds that are several times 10°M,, 
already in place before z = 14. GN-z11 would fall in this category, and 
the tracks obtained in ref. 3 would also explain the black hole to stellar 
mass ratio observed in GN-z11. 

Other studies have proposed other scenarios, visualizing different 
seeding mechanisms, at different redshifts, and with different assump- 
tions about the accretion and merging rates, and which are capable of 
reproducing the black hole mass of GN-z11 by z= 10.6, and generally 
also its black hole to stellar mass ratio*”°?7", 

In sum, the properties of the black hole in GN-z11 can be explained 
using different assumptions, which can be broadly grouped in heavy 
seeds accreting at sub-Eddington rates, or intermediate-light seeds 
experiencing super-Eddington phases and/or modelled with a boosted 
Bondi accretion. 

More statistics on objects such GN-z11 are required to discriminate 
between different scenarios. For the time being, to our knowledge, 
GN-z11 remains the most luminous object at z > 10 in all HST Deep fields 
(including CANDLES and Frontier Fields). It is hoped that JWST observa- 
tions on larger areas (for example, in Cosmos-WEB) will find more AGN 
at z>10 similar to GN-z11. For the time being, as discussed in the text, it 
is interesting to note that models and simulations were expecting a few 
accreting black holes with masses in the range 10°-10’M, at 10 <z < 11 
inthe JADES Medium-Deep survey in the GOODS fields*"®. Therefore, 


the discovery of a2 x 10°M, black hole in GN-z11 is not unexpected, 
and a few more might be present (probably accreting at a lower rate) 
inthe GOODS fields. 


GN-z11 and its large-scale environment 
We have shown that the high nitrogen enrichment of GN-z11is probably 
restricted to the BLR, whose small mass and compact size has probably 
undergone very rapid chemical enrichment, requiring only a few SNe. 
We note that the high chemical enrichment of GN-z11is not in contrast 
with the recent claim of pristine gas in the halo of GN-z11“. These claims 
are on totally different scales, with the pristine gas found several kpc 
away from GN-z11, whereas the high chemical enrichment is estimated 
to bein the nucleus of GN-z11. Regarding the claim of pristine gas inthe 
halo of GN-z11, models expect that high-z massive galaxies may host 
pockets of pristine gas in their haloes, even down toz = 3 (refs. 106,107). 


Data availability 


The electronic version of the processed data used to produce the 
figures (including the 1D and 2D spectra of GN-z11) is available at the 
JADAES web site (jades-survey.github.io/). The NIRSpec raw data can 
be accessed at the JWST archive (archive.stsci.edu). 


38. Eisenstein, D. J. et al. Overview of the JWST Advanced Deep Extragalactic Survey (JADES). 
Preprint at arxiv.org/abs/2306.02465 (2023). 

39. Jakobsen, P. et al. The Near-Infrared Spectrograph (NIRSpec) on the James Webb Space 
Telescope. |. Overview of the instrument and its capabilities. Astron. Astrophys. 661, A80 
(2022). 

40. Ferruit, P. et al. The Near-Infrared Spectrograph (NIRSpec) on the James Webb Space 
Telescope. II. Multi-object spectroscopy (MOS). Astron. Astrophys. 661, A81 (2022). 

41. Böker, T. et al. In-orbit performance of the near-infrared spectrograph NIRSpec on the 
James Webb Space Telescope. Publ. Astron. Soc. Pac. 135, 038001 (2023). 

42. Bushouse, H. et al. JWST Calibration Pipeline. Zenodo zenodo.org/records/10569856 
(2023). 

43. Planck Collaboration. Planck 2018 results: VI. Cosmological parameters. Astron. Astrophys. 
641, A6 (2020). 

44. Maiolino, R. et al. JWST-JADES. Possible Population III signatures at z = 10.6 in the halo of 
GN-z11. Preprint at arxiv.org/abs/2306.00953 (2023). 

45. Laha, S. et al. lonized outflows from active galactic nuclei as the essential elements of 
feedback. Nat. Astron. 5, 13-24 (2021). 

46. Rodriguez Hidalgo, P. High Velocity Outflows in Quasars. PhD thesis, Univ. Florida (2009). 

47. Hamann, F. & Sabra, B. The Diverse Nature of Intrinsic Absorbers in AGNs. In AGN Physics 
with the Sloan Digital Sky Survey (eds Richards, G. T. & Hall, P. B.) Vol. 311 of Astronomical 
Society of the Pacific Conference Series, 203 (Astronomical Society of the Pacific, 2004). 

48. Chen, Z.-F. & Pan, D.-S. Collective properties of quasar narrow associated absorption 
lines. Astrophys. J. 848, 79 (2017). 

49. Byun, D., Arav, N. & Hall, P. B. The farthest quasar mini-broad absorption line outflow from 
its central source: Very Large Telescope/UVES Observation of SDSS JO242+0049. 
Astrophys. J. 927, 176 (2022). 

50. Pettini,M. et al. New observations of the interstellar medium in the Lyman Break Galaxy 
MS 1512-cB58. Astrophys. J. 569, 742-757 (2002). 

51. Bischetti, M. et al. Suppression of black-hole growth by strong outflows at redshifts 5.8-6.6. 
Nature 605, 244-247 (2022). 

52. Maiolino, R. et al. Extreme gas properties in the most distant quasars. Astron. Astrophys. 
420, 889-897 (2004). 

53. Diamond-Stanic, A. M. et al. High-velocity outflows without AGN feedback: Eddington- 
limited star formation in compact massive galaxies. Astrophys. J. Lett. 755, L26 (2012). 

54. Diamond-Stanic, A. M. et al. Compact Starburst Galaxies with FastOutflows: Central 
Escape Velocities and Stellar Mass Surface Densitiesfrom Multiband Hubble Space 
Telescope Imaging. Astrophys. J.912, 11(2021). 

55. Sell, P. H. et al. Massive compact galaxies with high-velocity outflows: morphological 
analysis and constraints on AGN activity. Mon. Not. R. Astron. Soc. 441, 3417-3443 
(2014). 

56. Perrotta, S. et al. Physical properties of massive compact starburst galaxies with extreme 
outflows. Astrophys. J. 923, 275 (2021). 

57. Arav, N. et al. Evidence that 50% of BALQSO outflows are situated at least 100 pc from the 
central source. Astrophys. J. 857, 60 (2018). 

58. Xu, X., Arav, N., Miller, T. & Benn, C. VLT/X-Shooter Survey of BAL quasars: large distance 
scale and AGN feedback. Astrophys. J. 876, 105 (2019). 

59. Walker, A., Arav, N. & Byun, D. High mass flow rate in a BAL outflow of quasar SDSS J1130 + 
0411. Mon. Not. R. Astron. Soc. 516, 3778-3785 (2022). 

60. Chisholm, J., Prochaska, J. X., Schaerer, D., Gazagnes, S. & Henry, A. Optically thin spatially 
resolved Mgll emission maps the escape of ionizing photons. Mon. Not. R. Astron. Soc. 
498, 2554-2574 (2020). 

61. Kuraszkiewicz, J. K. et al. Emission line properties of active galactic nuclei froma 
post-COSTAR Hubble Space Telescope Faint Object Spectrograph spectral atlas. 
Astrophys. J. Suppl. Ser. 150, 165-180 (2004). 

62. Nagao, T., Maiolino, R. & Marconi, A. Gas metallicity in the narrow-line regions of 
high-redshift active galactic nuclei. Astron. Astrophys. 447, 863-876 (2006). 


63. 


64. 


65. 


66. 


67. 


68. 


69. 


70. 


71. 


72. 


73. 


74. 


75. 


76. 


77. 


78. 


79. 


80. 


81. 


82. 


83. 


84. 


85. 


86. 


87. 


88. 


89. 


90. 


91. 


92. 


93. 


94. 


95. 


96. 


97. 


98. 


99. 


Cleri, N. J. et al. Using [Ne V]/[Ne III] to understand the nature of extreme-ionization 
galaxies. Astrophys. J. 953, 10 (2023). 

Vanzella, E. et al. The unusual N Iv]-emitter galaxy GDS JO33218.92-2753027: star 
formation or AGN-driven winds from a massive galaxy at z=5.56. Astron. Astrophys. 513, 
A20 (2010). 

Glikman, E., Djorgovski, S. G., Stern, D., Bogosavljević, M. & Mahabal, A. Discovery of two 
spectroscopically peculiar, low-luminosity quasars at z~ 4. Astrophys. J. 663, L73-L76 
(2007). 

Gutkin, J., Charlot, S. & Bruzual, G. Modelling the nebular emission from primeval to 
present-day star-forming galaxies. Mon. Not. R. Astron. Soc. 462, 1757-1774 (2016). 
Nakajima, K. & Maiolino, R. Diagnostics for Poplll galaxies and direct collapse black holes 
in the early universe. Mon. Not. R. Astron. Soc. 513, 5134-5147 (2022). 

Ferland, G. J. et al. The 2017 Release Cloudy. Rev. Mex. Astron. Astrofis. 53, 385-438 
(2017). 

Eldridge, J. J. et al. Binary Population and Spectral Synthesis Version 2.1: construction, 
observational verification, and new results. Publ. Astron. Soc. Aust. 34, e058 (2017). 
Salpeter, E. E. The luminosity function and stellar evolution. Astrophys. J. 121, 161-167 
(1955). 

Shakura, N. I. & Sunyaev, R. A. Black holes in binary systems. Observational appearance. 
Astron. Astrophys. 24, 337-355 (1973). 

Capellupo, D. M., Netzer, H., Lira, P., Trakhtenbrot, B. & Mejia-Restrepo, J. Active galactic 
nuclei at z~1.5 - |. Spectral energy distribution and accretion discs. Mon. Not. R. Astron. 
Soc. 446, 3427-3446 (2015). 

Leighly, K. M. & Moore, J. R. Hubble Space Telescope STIS ultraviolet spectral evidence of 
outflow in extreme narrow-line Seyfert 1 galaxies. |. Data and analysis. Astrophys. J. 611, 
107-124 (2004). 

Maiolino, R., Salvati, M., Marconi, A. & Antonucci, R. R. J. The Ly-edge paradox and the 
need for obscured QSOs. Astron. Astrophys. 375, 25-29 (2001). 

Maiolino,R., Juarez, Y., Mujica, R., Nagar, N. M. & Oliva, E. Early star formation traced by 
the highest redshift quasars. Astrophys. J. 596, L155-L158 (2003). 

de Rosa, G. et al. Evidence for non-evolving Fe !I/Mg || ratios in rapidly accreting z~6 QSOs. 
Astrophys. J. 739, 56 (2011). 

Mazzucchelli, C. et al. Physical properties of 15 quasars at z> 6.5. Astrophys. J. 849, 91 
(2017). 

Shin, J., Nagao, T., Woo, J.-H. & Le, H. A. N. The Fe 1I/Mg 1! flux ratio of low-luminosity 
quasars at z~3. Astrophys. J. 874, 22 (2019). 

Sameshima, H. et al. Mg Il and Fe II fluxes of luminous quasars at z~2.7 and the evaluation 
of the Baldwin effect in the flux-to-abundance conversion method for quasars. Astrophys. 
J. 904, 162 (2020). 

Maiolino, R. & Mannucci, F. De re metallica: the cosmic chemical evolution of galaxies. 
Astron. Astrophys. Rev. 27, 3 (2019). 

Oesch, P. A. et al. The most luminous z~ 9-10 galaxy candidates yet found: the luminosity 
function, cosmic star-formation rate, and the first mass density estimate at 500 Myr. 
Astrophys. J. 786, 108 (2014). 

Oesch, P. A. et al. A remarkably luminous galaxy at Z=11.1 measured with Hubble Space 
Telescope grism spectroscopy. Astrophys. J. 819, 129 (2016). 

Ai, Y. L. et al. A comparative study of optical/ultraviolet variability of narrow-line Seyfert 1 
and broad-line Seyfert 1 active galactic nuclei. Astron. J. 145, 90 (2013). 

Xue, Y. Q. et al. The 2 Ms Chandra Deep Field-North Survey and the 250 ks extended 
Chandra Deep Field-South Survey: Improved point-source catalogs. Astrophys. J. Suppl. 
Ser. 224, 15 (2016). 

Vasudevan, R. V. & Fabian, A. C. Piecing together the X-ray background: bolometric 
corrections for active galactic nuclei. Mon. Not. R. Astron. Soc. 381, 1235-1251 (2007). 
Buisson, D. J. K., Fabian, A. C. & Lohfink, A. M. Coronal temperatures of the AGN ESO 
103-035 and IGR 2124.7+5058 from NuSTAR observations. Mon. Not. R. Astron. Soc. 481, 
4419-4426 (2018). 

Netzer, H. Bolometric correction factors for active galactic nuclei. Mon. Not. R. Astron. 
Soc. 488, 5185-5191 (2019). 

Willott, C. J. et al. Eddington-limited accretion and the black hole mass function at 
redshift 6. Astron. J. 140, 546-560 (2010). 

Trakhtenbrot, B., Netzer, H., Lira, P. & Shemmer, O. Black hole mass and growth rate at 
Z~4.8: a short episode of fast growth followed by short duty cycle activity. Astrophys. J. 
730, 7 (2011). 

Shen, Y. et al. A catalog of quasar properties from Sloan Digital Sky Survey Data Release 7. 
Astrophys. J. Suppl. Ser. 194, 45 (2011). 

Shen, Y. et al. Gemini GNIRS near-infrared spectroscopy of 50 quasars at z>5.7. Astrophys. J. 
873, 35 (2019). 

Pensabene, A. et al. The ALMA view of the high-redshift relation between supermassive 
black holes and their host galaxies. Astron. Astrophys. 637, A84 (2020). 

Harikane, Y. et al. A JWST/NIRSpec first census of broad-line AGNs at z=4-7: Detection of 
10 faint AGNs with Mpu- 10°-10® M,un and their host galaxy properties. Astrophys. J. 959, 
39 (2023). 

Mezcua, M. et al. Overmassive black holes in dwarf galaxies out to z~0.9 in the VIPERS 
survey. Astrophys. J. Lett. 943, L5 (2023). 

Vestergaard, M. & Osmer, P. S. Mass functions of the active black holes in distant 
quasars from the Large Bright Quasar Survey, the Bright Quasar Survey, and the 
color-selected sample of the SDSS Fall Equatorial Stripe. Astrophys. J. 699, 800-816 
(2009). 

Kormendy, J. & Ho, L. C. Coevolution (or not) of supermassive black holes and host 
galaxies. Annu. Rev. Astron. Astrophys. 51, 511-653 (2013). 

Maiolino, R. et al. JADES. The diverse population of infant Black Holes at 4<z<11: merging, 
tiny, poor, but mighty. Preprint at arxiv.org/abs/2308.01230 (2023). 

Marconi, A. et al. The effect of radiation pressure on virial black hole mass estimates and 
the case of narrow-line Seyfert 1 galaxies. Astrophys. J. 678, 693-700 (2008). 

Khatu, V. C. et al. Supermassive black holes with high accretion rates in active galactic 
nuclei. XIII. Ultraviolet time lag of HB emission in Mrk 142. Astrophys. J. 958, 127 
(2023). 


Article 


00. Kocevski, D. D. et al. Hidden little monsters: Spectroscopic identification of low-mass, 
broad-line AGNs at z > 5 with CEERS. Astrophys. J. Lett. 954, L4 (2023). 

01. Ding, X. et al. Detection of stellar light from quasar host galaxies at redshifts above 6. 
Nature 621, 51-55 (2022). 

02. Scholtz, J. et al. GN-z11: the environment of an AGN at z=10.603. Preprint at arxiv.org/ 
abs/2306.09142 (2023). 

03. Zhang, H. et al. TRINITY I: self-consistently modelling the dark matter halo-galaxy- 
supermassive black hole connection from z=0-10. Mon. Not. R. Astron. Soc. 518, 
2123-2163 (2023). 

04. Smidt, J., Whalen, D. J., Johnson, J. L., Surace, M. & Li, H. Radiation hydrodynamical 
simulations of the first quasars. Astrophys. J. 865, 126 (2018). 

05. Trinca, A. et al. Seeking the growth of the first black hole seeds with JWST. Mon. Not. R. 
Astron. Soc. 519, 4753-4764 (2023). 

06. Liu, B. & Bromm, V. When did Population III star formation end? Mon. Not. R. Astron. Soc. 
497, 2839-2854 (2020). 

07. Venditti, A. et al. A needle in a haystack? Catching Population III stars in the epoch of 
reionization: |. Population III star-forming environments. Mon. Not. R. Astron. Soc. 522, 
3809-3830 (2023). 


Acknowledgements We acknowledge the suggestions and discussions we had with several 
colleagues, in particular: G. Risaliti, A. Marconi, J. Bennett, S. Koudmani, D. Sijacki, R. Schneider, 
M. Volonteri, R. Valiante, A. Trinca, D. Berg, R. Ellis, J. Dunlop, M. Pettini, H. Katz, W. N. Brandt 
and R. Tripodi. F.D., J.S., L.S., R.M. and T.J.L. acknowledge support from the Science and 
Technology Facilities Council (STFC), from the ERC through advanced grant 695671 ‘QUENCH’, 
and from the UKRI Frontier Research grant RISEandFALL. R.M. also acknowledges funding from 
a research professorship from the Royal Society. A.B., G.C.J. and J.C. acknowledge funding 
from the ‘FirstGalaxies’ Advanced Grant from the European Research Council (ERC) under the 
Horizon 2020 research and innovation programme of the European Union (grant agreement 
No. 789056). B.R. acknowledges support from the NIRCam Science Team contract to the 
University of Arizona, NAS5-02015. B.R.D.P., M.P. and S.A. acknowledge support from grant 
PID2021-127718NB-I00 funded by the Spanish Ministry of Science and Innovation/State Agency 
of Research (MICIN/AEI/10.13039/501100011033). M.P. also acknowledges support from the 
Programa Atraccion de Talento de la Comunidad de Madrid by grant 2018-T2/TIC-11715. 
C.N.A.W., E.E. and F.S. acknowledge a JWST/NIRCam contract to the University of Arizona 
NAS5-02015. D.J.E. is supported as a Simons Investigator and by JWST/NIRCam contract to the 


University of Arizona, NAS5-02015. E.C.-L. acknowledges support of an STFC Webb Fellowship 
(ST/W001438/1). H.Ü. acknowledges support from the Isaac Newton Trust and from the Kavli 
Foundation through a Newton-Kavli Junior Fellowship. J.W. acknowledges support from the 
ERC advanced grant 695671, ‘QUENCH’, and the Foundation MERAC. S.Carniani acknowledges 
support by HE ERC Starting grant no. 101040227 - WINGS of the European Union. The research 
of C.C.W. is supported by NOIRLab, which is managed by the Association of Universities for 
Research in Astronomy (AURA) under a cooperative agreement with the National Science 
Foundation. This research is supported in part by the Australian Research Council Centre of 
Excellence for All Sky Astrophysics in 3 Dimensions (ASTRO 3D), through project no. 
CE170100013. 


Author contributions R.M., J.S., JW., S.C., F.D., A.d.G., H.Ü. and ST. contributed to the writing of 
the paper, methods and creation of figures. All authors contributed to the interpretation of the 
results. N.K., B.R.d.P. contributed to the design, construction and commissioning of NIRSpec. 
S.A., S. Carniani, M.C., J.W., M.P. and B.R.d.P. contributed to the NIRSpec data reduction and 

to the development of the NIRSpec pipeline. S.A. contributed to the design and optimization 
of the MSA configurations. A.B., C.N.AW., E.C.-L., K.B. and H.U. contributed to the selection, 
prioritization and visual inspection of the targets. S. Charlot, J.C., E.C.-L., R.M., J.W., F.D., T.J.L., 
M.C., A.d.G. and L.S. contributed to the analysis of the spectroscopic data, including redshift 
determination and spectral modelling. F.D., T.J.L., M.C., B.R.d.P., R.M., S.A. and J.S. contributed 
to the development of the tools for the spectroscopic data analysis, visualization and fitting. 
CW. contributed to the design of the spectroscopic observations and MSA configurations. 
C.N.A.W., C.J.W., D.J.E., R.M. and S.A. contributed to the design of the JADES survey. E.E., K.N.H. 
and C.CW. contributed to the design, construction and commissioning of NIRCam. B.R., D.J.E., 
I.S., ST., C.N.A.W. and Z.J. contributed to the JADES imaging data reduction. B.R. contributed to 
the JADES imaging data visualization. 


Competing interests The authors declare no competing interests. 


Additional information 

Supplementary information is available for this paper at https://doi.org/10.1038/s41586-024- 
07052-5. 

Correspondence and requests for materials should be addressed to Roberto Maiolino. 

Peer review information Nature thanks the anonymous reviewers for their contribution to the 
peer review of this work. 

Reprints and permissions information is available at http://www.nature.com/reprints. 


observed wavelength (um) 


MEN 3.20 3.22 3.24 3.26 3.28 
k 
jy 
[= 
o 
T 
a 
D 
o 
o 
b 
< 
a 2,760 2,780 2,800 2,820 2,840 
rest-frame wavelength (A) 
observed wavelength (um) 
= 1.41 142 143 144 145 1.46 
k 
T 
E 
o 
T 
oo 
D 
o 
3 
o 
< 
a 1,220 1,230 1,240 1,250 1,260 
rest-frame wavelength (A) 
o O20by BE i 3 
N Sh to i2 
g 0.15 = i22 i= 
> z 12 © 10 
4 © N 85 iz 
o 0.10 Z 13 © i~ 
D m= 
o 0.05 LA 
a i . . 
a 0.00 Refer kd fase 
i ooh 


3,850 3,875 3,900 3,925 3,950 3,975 


observed wavelength (um) 


—< 1.88 1.90 1.92 1.94 
< 7 
y {Hell (1640) 
§ 
T 
o 
D 
o 
o 
T 
= 
z 
a 1,620 1,640 1,660 
rest-frame wavelength (A) 
observed wavelength (um) 
æ 218 2.20 2.22 2.24 
k 
* o2 Ed) cim (1906) 'x CIN] (1908) 
§ 
1 
9 01 
o 
S 
a d j RED Ea E ed dl da aadi ELER -a LEETETTT 
2 0.0 
= ! . 
1,880 1,900 1,920 
rest-frame wavelength (A) 
i 
i 
| Hô 
I 
1 


4,075 4,100 4,125 4,325 4,350 4,375 


rest-frame wavelength (A) 


Extended Data Fig. 1| Zoomin onthe additional emission lines fitted. 

a) MgIIA2796,2804 doublet; b) HellA1640; c) Lya, NVA1238,1242 doublet 
(undetected) and SilIA1260,1264 (undetected), corrected for the Lya damping 
wing; d) C111]A1906,1908 doublet. As the doublet is unresolved, the fit turns out 
degenerate between line width and fluxes of the two components; moreover it 


is also contributed to by star formation in the host galaxy (see text for details); 
e) [Ne111]A3869 profile compared with the Balmer lines Hô and Hy. In all panels, 
the continuum is subtracted. The black dotted lines indicate the lsigma noise 
level. 
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Extended Data Fig. 2 | NV/CIV versus NV/Hell flux ratio diagram. GN-z11 
(golden circle) is compared with the ratios observed for the broad lines of type 
1AGN (purple stars”), and for the NLR of type 2 AGN (red squares”), illustrating 
that the non-detection of NV for GN-z11is not constraining and consistent with 
the AGN scenario. 
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Extended Data Fig. 3 | CIII]/CIV versus CIII]/Hell flux ratio diagrams. GN-z11 
(golden circle) is compared with: a) the AGN-NLR models (red squares) by? and 
SF galaxies (blue stars) by® (left) and with b) the AGN-NLR models (red squares) 
and SF models (blue triangles) by” (centre). All models have been chosen ina 
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low metallicity range, around the value inferred by’ for GN-z11. c) Comparison 
of GN-z11 with the ratios observed for the broad lines of type 1 AGN (purple 
circles”), narrow lines of type 2 AGN (red squares”), and starburst galaxies 
(blue stars'’), 
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Extended Data Fig. 4 | Flux ratios of density-sensitive nitrogen lines as a function of hydrogen gas density, n,,. Same as Fig. 2 but where we have separated the 
photoionization models for AGN (left) and Star Forming galaxies (right). 
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Extended Data Fig. 5| Low resolution (prism) spectrum of GN-z11. Left: 
Observed prism spectrum (black solid) compared with the (maximum) 
contribution from the host galaxy of the AGN as inferred by’ (orange dotted), 
and the nebular emission inferred from a simple Cloudy model (purple dashed) 
normalized to the Hy flux not included in the galaxy model. Right: Spectrum 
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subtracted of the galactic and nebular continua, ina log-log scale, whose regions 
not affected by emission lines (solid black) have been fitted with a simple 
powerlaw (red-dashed) resulting into a slope of -2.26 + 0.10, consistent with 
the slope expected for an accretion disc (—2.33, dotted green line). 
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Extended Data Table 1| Measured fluxes and FWHM of the emission lines used in this work 


Emission line Ao Flux Doublet-Multiplet FWHM 
A (x10719) ergs s~! cm~? ratios km s~! 
NV 1239+1243 <3.50 
[NIV] 1483.32 < 1.6% 
CII* 1335 5627) 380 + 80 
NIV] 1486.50 10.4245 470 + 50 
CIV (emission) 1550.77 SPAA 
Hell 1640.50 6:00 355 1200 + 380 
NII] multiplet-total 10:37 45 430 + 80° 
1746.82° /total 0.025705 
1748.64¢ /total Goan ce 
1749.67 /total 0:3700. 
1752.16°/total 0.14+0:06 
1753.99% /total 02528 
CuI] 1906+1908 9.9775 610 + 60° 
Mell 2796+2804 46008 430 + 65/ 
2796/2804 1s ear 
[NeV] 3426.86 <1.5¢ 
[Nerv] 2424+2422 (fit)9 7 a 380 + 110 
242442422 (int)9 SI 
[o1] 37274+3729 ONEEN 368 + 50 
3729/3727 0:59:57 
[NelI]] 3869.68 LOE oe 340 + 30 
Hel+H8 3889.73 State 360 + 50 
Hô 4101.73 Se 410 + 70 
Hy 4340.47 yey 360 + 80 


Notes: °30 upper limit (in the case of NV is the upper limit on the sum of the doublet); *the four components of the multiplet are forced to have the same width in the fit; ĉas the components 
M1746.82 and A1752.16 come from the same level, their fluxes are forced to have the ratio given by their Einstein's coefficients, i.e. 1746.82/1752.16=0.146; das the components A1748.64 and 
41753.99 come from the same level, their fluxes are forced to have the ratio given by their Einstein's coefficients, i.e. 1753.99/1748.64=1.05; “for the CIII] doublet is totally blended and is here 
fitted with a single Gaussian (see text for a more complex, but uncertain decomposition); ‘the two components of the Mgll doublet are forced to have the same width; Sin the case of the [NelV] 
doublet, as it is a key line, we provide both the result from MCMC fitting, assuming a single Gaussian (first line) and also by simply taking the integrated flux (second line). 


Extended Data Table 2 | Features observed in GN-z11 versus possible energy sources 


Source 
AGN W olf-Rayets N uclear 
Feature (NLSy1) Starburst 
(without WR) 

Gas density > 10° cm’ 3 Consistent Inconsistent Inconsistent 
(from NIII & NIV lines) 

NelV detection Consistent Inconsistent Inconsistent 
CII’ detection Consistent Inconsistent Inconsistent 
Hell broad Consistent Consistent Inconsistent 
NelV (1718)/ Nel V(1485)< 10| Consistent Inconsistent Consistent 
Vwing > 800km/s Consistent Inconsistent Inconsistent 
UV-opt. cont. slope = 2.3 Consistent Consistent Consistent 
X-ray upper limit Consistent Consistent Consistent 
CIII]/ CIV vs CIII]/ Hell Consistent Consistent Consistent 
NV upper limit Consistent Consistent Consistent 


